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In this paper, we study axially symmetric solutions with i? = 2 — 5 in the chiral quark soliton 
model. In the background of axially symmetric chiral fields, the quark eigenstates and profile 
functions of the chiral fields are computed self-consistently. The resultant quark bound spectrum 
are doubly degenerate due to the symmetry of the chiral field. Upon quantization, various observable 
spectra of the chiral solitons are obtained. Taking account of the Finkelstein- Rubinstein constraints, 
we show that our results exactly coincide with the physical observations for B = 2 and 4 while B = 3 
and 5 do not. 
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I. INTRODUCTION 

The chiral quark soliton model(CQSM) was developed 
in 1980's as a low-energy effective theory of QCD. Since 
it includes the Dirac sea quark contribution and ex- 
plicit valence quark degrees of freedom, the model in- 
terpolates between the constituent quark model and the 
Skyrme model 0,11111113. The CQSM is derived from 
the instanton liquid model of QCD vacuum and incor- 
porates the non-perturbative feature of the low-energy 
QCD, spontaneous chiral symmetry breaking. It has 
been shown that the B — 1 solution provides correct 
observables as a nucleon including mass, electromagnetic 
value, spin carried by quarks, parton distributions and 
octet SU(3) baryon spectra. 

For B = 2, the stable axially symmetric soliton solu- 
tion was found in |^. The solution exhibits doubly de- 
generate bound spectrum of the quark orbits in the back- 
ground of the axially symmetric chiral field with winding 
number two. Upon quantization, various dibaryon spec- 
tra were obtained, showing that the quantum number of 
the ground state exactly coincide with that of physical 
deutcron nil. For B > 2, the Skyr me model predicts 
that the solutions have only discrete, crystal-like symme- 
tries H 0, 0| . According to the prediction, we stud- 
ied the CQSM with B = 3 tetrahedrally symmetric chi- 
ral fields and obtained triply degenerate spectrum of the 
quark orbits . Its large degeneracy indicates that the 
terahedrally symmetric solution may be the lowest-lying 
configuration. For B > 3, one can also expect that the 
lowest-lying solutions in the CQSM inherits the discrete 
symmetries predicted in the Skyrme model. Studying 
solutions with those symmetries in CQSM is, however, 
formidable especially for quantization. Thus, before em- 
barking those discrete symmetries, it will be instructive 
to study axially symmetric solutions which are much sim- 



pler. Besides, considering the fact that for some higher 
baryon numbers, the ground states of the skyrmions do 
not agree with the experimental observation jl3] , the 
possibility that axially symmetric solutions provide cor- 
rect ground states can not be excluded. Recently, it was 
found in Ref. that in the BPS monopoles all axially 
symmetric solutions up to chrage five have lower energy 
than that of discrete symmetries. A cylindrical shape 
isomer in ^^C was also found in the Skyrme model frame- 
work . Research of axially symmetric solitons is thus 
in progress. 

In Sec^ we shall obtain axially symmetric classical 
soliton solutions with B = 2 — 5. The solutions exhibit 
doubly degenerate spectra due to their axial symmetry. 
Such degeneracy generates large shell gaps and confirms 
that the solutions are stable local minima. In Sec lIIII we 
shall quantize the obtained classical solitons semiclassi- 
cally. Imposing the Finkelstein-Rubinstein constraints 
on the states, the ground states of the axially symmetric 
solitons are constructed and examined if they agree with 
the experimental observation. In Sec II VI is the detail of 
the numerical analysis used to obtain the classical and 
quantum solutions. Conclusions and discussions arc in 
SeclVl 



II. AXIALLY SYMMETRIC CLASSICAL 
SOLUTIONS 

The CQSM is derived from the instanton liquid model 
of the QCD vacuum and incorporates the nonperturba- 
tive feature of the low-energy QCD, spontaneous chiral 
symmetry breaking. The vacuum functional is defined 

by i 



Vir'Dtp'Dip^ exp 



where the SU(2) matrix 
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describes chiral fields, ■tp is quark fields and M is the 
constituent quark mass, f^^ is the pion decay constant 
and experimentally /^r ~ 93MeV. 

The B — 1 soliton solution has been studied in de- 
tail at classical and quantum level in P, 0, 0, 0, ■ To 
obtain solutions with _B > 1, we shall employ the chi- 
ral fields with winding number B in the Skyrme Model 
as the background of quarks, which can be justified as 
follows. 

In Eq.(^), performing the functional integral over tp 
and V'^ fields, one obtains the effective action 

S'eff(C/) = -iNcSplniD = -iNJogdetiD, (2) 

where iD — MU^'-" is the Dirac operator. The classi- 
cal solutions can be obtained by the extremum condition 
of l(2Jl with respect to U . For this purpose, let us consider 
the derivative expansion of the action 0, 0| . Up to 
quartic terms, we have, 



where 



N. 



CtriL^L^ 



-^tr{l[L„L.]^4(a,L^) + l(L,L^f} 



(3) 



where — d^UW . Suitably adjusting the coeffi- 
cients, one can identify ^ with the Skyrme model action. 
Therefore, it will be justified to adopt the configurations 
of the solutions in the Skyrme model to chiral fields in 
the CQSM. 

In the Skyrme model the minimal energy pion config- 
uration with B — 2 has an axial symmetry |l8| and can 
be written by 



U{x) — cos F{p, z) + iT ■ fi sin F{p, z) 



(4) 



where 



n — (sin z) cosm^Lp, sin z) sinm„(^, cos 6(/5, z)) 

(5) 

and mw is the winding number of the pion fields. We 
shall use this configuration in the backgound to obtain 
axially symmetric chiral quark solitons. 

In the CQSM, the number of valence quark is associ- 
ated with the baryon number such that the baryon num- 
ber B soliton consist oi NcX B valence quarks. If the cor- 
relation between quarks is sufficiently strong, their bind- 
ing energy become large and the valence quarks can not 
be observed as positive energy particles |23,|2J|. Thus, 
one gets the picture of the topological soliton model in the 
sense that the baryon number coincide with the winding 
number of the background chiral field when the valence 
quarks occupy all the levels diving into negative energy 
region. 

Let us rewrite the effective action in ^ as 



= -iNc\ogdet{idt - HiU^"-)) 



(6) 



(7) 



The classical energy of the soliton can be estimated from 
the quark determinant in Eq. © |2^ . We introduce 
the eigenstates of operates, idt — H{U^^) and H{U'^^), 
such that 



(8) 
(9) 



where 5'^ „ = e 



and A 



-E^ + ujn- Impos- 



ing on ^'^.n the anti-periodicity condition, \E'^.„(a;,T) 
-^'M,n(a;,0), reads 

w„T = (2n -f l)7r. (10) 

The determinant in Eq.© then becomes 

detiidt-H) = []A^,„ 



(2n + l)7r 



c n (i ^"'^'"^ 



/^,n>0 



(2n+l)27r2 



= cncos(i|i?jr) 
= f exp(zi^|i?,|r) 

xY[(l + exp{~i\E^\T)) (11) 



where 



(2n+l) V2 

y2 



and the product formula for the cosine function cos(z) — 
n^>i(l-4z2/(2n-l)27r2) has been used. Inserting l(TT|) 
into ©, one obtains 

5eff = -iV,r^n^|S^| + Ar,rJ^|£;^| (12) 



where is the valence quark occupation number which 
takes values only or 1. Correspondingly, the classical 
energy is given by 



E — E,,a] + El 



val -t- -C/fiold 



(13) 



where 



-Enfield — ^9^c^ \Ei 
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representing the valence quark and sea quark contribu- 
tion to the total energy respectively. 

The effective action SesiU) is ultraviolet divergent and 
hence must be re gula rized. Using the proper-time regu- 
larization scheme l27l. we can write 



where 



srm = r -sp (e-^'^- 



-D'Dor 



-NcT / — 



1/A2 T 



xSp 



(14) 



where Dq and Hq are operators with U — 1. The total 
energy is then given by 

Est.tic[U] = E,,,[U] + Afield [C/] - Afield [[/ = 1] (15) 

where 



val 



with 



A 



exp 



E^ 



\2' V A 



and E'^!^'} is the valence energy of the i th valence quark. A 
is a cutoff parameter evaluated by the condition that the 
derivative expansion of p4l) reproduces the pion kinetic 
term with the correct coefficient, i.e., 



f 



47r2 



dr 



The extremum conditions for the total energy 

5 ^ .... . S 



(16) 



SFip, z) 



E.. 



-Bstatic[C/] =0 (17) 



yield the following equations of motion for the profile 
functions, 

R^{p,z)cosQ{p,z) = R^{p,z)sme{p,z) (18) 
S{p, z) smF{p, z) = P{p, z) cos F{p, z) (19) 

where 

P{p,z) = i?^(/9, z)sine(p, z) -f i?^(p, z)cose(p, z) . 

(20) 

In terms of eigenfunction (f) in Eq.ljS]), and S are 

given by 

R^{p,z)^R^^,ip,z) + R^{p,z) (21) 
i?^(p,z) = i?ii(p,z)-hi?o^(p,z) (22) 
S{p,z) = Svi,i{p,z) + So{p,z) (23) 



i 

X (ti cosm^ip + T2 smm^ip)(f>i{p, (p, z) , 
RoiPiz) = ^J^iEf,)sgn{E^) / dip(j)f,{p,(p,z)i"f5 

x(ti cos rriyjip + T2 sinmwV')'/'M(P' ^) ' 

^val(P'^)=X! / ^'/'^«(P''/^'^)*75T3(/'»(p,(p,z) , 

RoiP^z) = J2-^{Ef,)sgn{Ef,) / dip4>f,{p,tp,z)i-/5 

'S'val(p,z) = ^ / dip4),{p,^p,z)(l)i{p,ip,z) , 
i 

So{p,z) ^^AfiE^)sgn{Ef,) J dipct)f,{p,ip, z)ct)f,{p,ip, z) . 



and subscripts, and val, represent the vacuum and va- 
lence quark contributions respectively. The boundary 
conditions for the profile functions were constructed by 
Braaten and Carson ; 

F{p,z)^0 as p^ + z^— >cx), 

i^(0,0) = -7r, e(0,z)=|||'^ l>l. (24) 

The procedure to obtain the self-consistent solutions of 
Ea. lfH^ and (fH^ is 1) solve the eigenequation in Eq.Q 
under assumed initial profile functions Fo(p, z), 0o(p, z) 
which satisfy the boundary conditions eas. (|24|l . 2) use 
the resultant eigenfunctions and eigenvalues to calculate 
R^,R^,S and also P in Eq.^, 3) solve Eq. JTHl and 
(|19|1 to obtain new profile functions, and 4) repeat 1) — 3) 
until the self-consistency is attained. 

In Figs. I1I4I we show the spectrum of the quark orbits 
in the background of chiral fields with winding number 
Tow = 2 — 5, as a function of the size parameter X. The 
profile functions are parameterized by X as 



F{p,z) = -TT + n^/p^ + z'^/X for X < ^/p^~^ 

= otherwise, (25) 

e(/9,z) = tan"i(p/z). (26) 

To examine the spectrum in detail, let us consider the 
Hamiltonian defined in O. For the axially symmetric 
chiral field in Eq.Q, this Hamiltonian commutes with 
the third component of the grand spin operator and 
the time-reversal operator T. These are specifically, 



1 



K3 = L3 

T = i7i73 • iriTaC 



1 



-0-3 + -TO„T3 



(27) 
(28) 
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FIG. 1: Spectrum of the quark orbits with B = 2 as a func- 
tion of the soliton size parameter X. 



where L3, 173, and T3 are respectively the third component 
of orbital angular momentum, spin, and isospin operator, 
and C is a charge conjugation operator. The parity op- 
erator is defined by P = 70 for odd -B, and P = 70 T3 
for even B. Thus, the eigenvalues of the Hamiltonian 
can be specified by the magnitude of K3 and the parity 
TT = ±. We have K3 = 0, ±1, ±2, ±3, • • • for odd B, and 
if3 = ±i, ±|, ±|, ±|, • • • for even B. Since the Hamil- 
tonian is invariant under time reverse, the states of +Kt, 
and —K3 are degenerate in energy. 

Figffl shows the quark spectrum with = 2. The 
bound states diving into negative region are doubly de- 
generate with K3 — ±i. Thus putting Nc — 3 va- 
lence quarks on each of the bound levels, we have the 
B = 2 soliton solution. For = 3, the spectrum 
of if 3 — ±1" (double degeneracy) and K3 — 0+ states 
dive into negative-energy region. Thus, we have the 
B — 3 soliton solution. For rn„ — 4, the spectrum 
oi K3 = and K3 = ±^ (both doubly degener- 

ate) states dive into negative region. Thus, we have the 
_B = 4 soliton solution. For ~ 5, the spectrum of 
Ks = ±2+(double), K3 = ±1" (double) and K3 = 0+ 
states dive into negative-energy region. Thus we have the 
B — 5 soliton solution. These results confirm that the 
baryon number of the soliton is identified with the num- 
ber of diving levels occupied by Nc valence quarks. It can 
be seen that the degeneracy which occurs due to symme- 
try of the chiral field reduces the number of states and 
hence makes large shell gaps. This observation indicates 
that degeneracy in spectrum contribute to make classical 
energies of the soliton solutions lower. In fact, our B = 2 
solution which is considered to be the minimum energy 
soliton from the study of the B = 2 skyrmions provides 
the maximum degeneracy in spectrum. It will be in- 
teresting to study minimum solutions from this point of 
view. 

The baryon number density is defined by the zeroth 
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FIG. 2: Spectrum of the quark orbits with B = 3. 
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FIG. 3; Spectrum of the quark orbits with B = 4. 
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FIG. 4: Spectrum of the quark orbits with B = 5. 
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FIG. 5: Contour plot of the profile functions with B = 2. 
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FIG. 6: Contour plot of the profile functions with B = 3. 
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FIG. 7: Contour plot of the profile functions with B = 4. 
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FIG. 8: Contour plot of the profile functions with B = 5. 
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TABLE I: The classical mass for B = 1 ~ 5 {[MeV]). 



where 



B 




Valence 


Vacuum Total 


1 173 






674 1192 


2 173 


173 




1166 2204 


3 173 


173 


298 


1561 3493 


4 106 


106 


232 232 


2727 4753 


5 145 


145 


319 319 


409 2537 6543 


TABLE IL 


The mean radius and mean root square radius in 


each baryon number for B — 1 ~ 


5 .(The result of B = 1 is 


quoted by Ref. 5] 


)• 




B 




{P)[fm] 




1 






0.785 


2 




0.672 


0.821 


3 




0.659 


0.854 


4 




0.971 


1.140 


5 




1.048 


1.225 



component of the baryon current 
1 



where 



(B) = — (^7oV^) = (B>.ai + (5)o (29) 



i 



5]AA(4°))sgn(4°)) 



The contour plot of the baryon number density for each 
baryon number is shown in Fig. |51 It can be seen that 
they have toroidal in shape. 

The mean radius (p) is given by 



(P) = (P)val + (P)0- 



(30) 



where 



(P>vai = — X! / Pdpdzd'fP4>\{P^V,z)(t>i{p,'P,z) 
X / pdpdzdLpp(j)Up,if,z)(j)f,{p,if,z). 



The root mean square radius is given by 

+ (31) 



(r^)vai = — X! / pdpdzdip{p^ + z'^)(f>l{p, ip, z)4,i{p, Lp, z) 



X / pdpdzdip{p'^ + z'^)(j)l{p,ip,z)(j)i,{p,ip,z) . 



These values for each baryon number are shown in Ta- 
ble llll 



III. QUANTIZATION 

A. Rotational Zero Mode Quantization 

The sohtons that we obtained in the previous section 
are classical objects and therefore must be quantized to 
assign definite spin and isospin to them. Quantization 
of the solitons can be performed semiclassically for their 
rotational zero modes. For the hedgehog soliton, because 
of its topological structure, a rotation in isospin space 
is followed by a simultaneous spatial rotation. For the 
axially symmetric soliton, there are six rotational zero 
modes by rotations of iso-degrees of freedom and spatial 
rotations. 

Let us introduce the "dynamically rotated" chiral fields 

M 

Uix,t) = A{t)Uix')A{t)\ x''^mB{t)]x^ (32) 
where 

E][B{t)]^lTr[a^Bit)a,B{t)^], (33) 

and A{t) and B{t) are time-dependent SU{2) matrices 
generating an iso-rotation and a spatial rotation respec- 
tively. By transforming the rotating frame of reference, 
the Dirac operator with Eq. H32|l can be written as 

i~D ^ i^- MU^'^{x,t) 

= A{t)Sit)^-/"[idt + i-^^'^'^dk' - V^x') + iX^A 
^iS^ S]S{t)A{t)^ (34) 

where 

^,0 



(35) 



and S{t) is the rotation operator for the Dirac field and 6 
is an angle of the spatial rotation. Note that the gamma 
matrices 7^ explicitly depend on the coordinates and do 
not transform as a contravariant vector p^ . Substitut- 
ing Eq. H35|l into Eq. (|34|l . one obtains 



i~D = A{t)S{t)^"f°[idt - H{W) + nA + nB]S{t)A{t)'^ 



(36) 
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FIG. 9: baryon number density [fm '^]. 



where 

nA = zAU = i^i^Ta (37) 
ng iS^S +{rxp)-e^ n%Jb (38) 

with Ja = l/2eabcl^^^ — «(»" X V)a- and Q.b are the 
angular velocity operators for an isorotation and for a 
spatial rotation respectively. Assuming that the rotation 
of the soliton is adiabatic, we shall expand the effective 
action S'e// around the classical solution Us with respect 
to the angular momentum velocity ^Ia and up to 
second order 

= S,s{Us) - iN, Sp [\og{idt - H{Uf) + ^A + ^b)] 
-S^[\og{idt-H{Uf)] (39) 



With the proper-time regularization, we have 

+ \ J di[/o^df^^(t)f^^(t)+/o^f,f7^(t)f7^B(i) 

+i^^^Mt)n\(t) + i^f,n%mUt)] (40) 

where IqS are the vacuum sea contributions to the mo- 
ments of inertia defined by 



tAA 
^O.ab — 


71, m 


Ta 


m) {'m\T}j\n) 


tAB 
^Q,ab — 


iiV,^/(i?™,i?„,A)(n 

n,m 


Ta 


m) {m\Jb\n) 


jBA 


iiV,^/(i?™,£;„,A)(n 


Ja 


m) (m|Tb|n) 


tBB 
^Q,ab — 


iiV,^/(i?™,i?„,A)(n 


Ja 


m)(TO Jf, n) 



7i,m 
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with the cutoff function f{Em., En, A) 

fiE„.,En,A) 

2A e-^^/^^ - e-^n/A^ 
- ^ El -El 

, sgn(£;„)erfc(|£;„|/A) - sgn(S„)ei-fc(|S„|/A) 



Er, 



(41) 



Similarly, for the valence quark contribution to the mo- 
ments of inertia, we have 



jAA 
val,a6 



jAB 
val,a6 



tBA 
val.ab 



tBB 
val,ab 



(val I Ta I m) (mirfel val) 



Eval 



(val|Ta|m)(m|Jj,| val) 
(val I Ja|m) (to I Tb I val) 
(val I Ja I m) (to I Jb I val) 



Evcd 



(42) 



The total moments of inertia are then given by the sum 
of the vacuum and valence as 



tAA 
^ab 


- T^"^ A 
val,a& ^ 


tAA 
~ ^0,ab 


tAB 
^ab 


— r^s J 

val,ab 


, tAB 
r ^0,ab 


tBA 
^ab 


— jBA 
val,a6 


, jBA 


tBB 
^ab 


— jBB 
val.afc 


, rBB 



From axial symmetry of the system, following relations 
are derived 



TABLE III: Moments of inertia. 



B 




Valence 


Sea 


Total 


2 


tAA 

-'ii 

rBB 

^11 

tAA 
33 


0.00773 
0.01141 
0.00429 


0.00363 
0.00464 
0.00125 


0.01136 
0.01605 
0.00554 


3 


rAA 

Jll 

rBB 
^11 

rAA 
J33 


0.01231 
0.02174 
0.00594 


0.00280 
0.00384 
0.00027 


0.01511 
0.02558 
0.00622 


4 


tAA 
-'11 

rBB 
-'11 

7-AA 
J 33 


0.01408 
0.04272 
0.01172 


0.00959 
0.01245 
0.00074 


0.02366 
0.05517 
0.01246 


5 


rAA 
-'11 

rBB 
^11 

rAA 
J 33 


0.02786 
0.12124 
0.01368 


0.00716 
0.01112 
0.00007 


0.03502 
0.13236 
0.01375 



The quantization conditions for the collective coordi- 
nates, A{t) and B{t), define a body-fixed isospin operator 
K and a body-fixed angular momentum operator L as 



rAAQb I rABQb 
^ab "a + -'ab 



jBAQb I jBBQb 
^ab "a + ^ab "S 



tr M 



2 dA 



-Ka (45) 



These are related to the usual coordinate-fixed isospin 
operator la and coordinate-fixed angular momentum 
operator by transformations. 



^6 r 



(47) 



To estimate the quantum energy corrections, let us in- 
troduce the basis functions of the spin and isospin oper- 
ators which were inspired from the cranking method for 
nuclei IHIp, 



jAA 



-'11 — ^22 



ij =0, 3, 

jAA jBB _ jBB 



-'11 



AB 



tAB 



^11 



jBB _ 2rAA jAB 
-'33 ~ '"■ -'33 ' ^33 



-'11 — ^22 1 
BA _ jBA 
^22 

BA _ 



'33 



0, 

TO„4^. (43) 



Theoretically, these moments of inertia can be computed 
using the eigenstates of Eq.(|5J). However, due to the dif- 
ference of the boundary conditions between the initial 
and final states of the matrix element, the moments of 
inertia acquire nonzero values with vanishing pion fields. 
To overcome this problem, we make the following replace- 
ment 13(1 



{n\Ja\m) 



{n\[H{U}'),Ja]\m)/{En-En 



En 



(44) 



where la — —i{rxV)a- Unless the Hamiltonian explicitly 
depend on the coordinates, the numerator vanishes with 
vanishing pion fields. The spurious contributions to the 
moment of inertia can be removed in this way. 



{A,B\ii3k3,jj3l3) 



V(2» + l)(2j + l) 
87r2 



where D is the Wigner rotation matrix. Then, we find 
the quantized energies of the soliton as 



E 



E. 



1 



OjAA 



1) 



1 



OjBB 



tAA 
^33 



tAA tBB 
^11 ^11 



k1 



(48) 



where i{i + + 1) and are eigenvalues of the 

Casimir operators and J^, and the operator K3, re- 
spectively. 

On Table Iml are the results of our calculation of mo- 
ments of inertia, and /^^, with B = 2-5. 
It is instructive to compare our results with the Skyrme 
model 22] where Un = 0.104, Vn = 0.163 and U33 = 
0.0709 which are correspondingly our Iii^ and I^/^- 
They are qualitatively in good agreement. 
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B. Finkelstein- Rubinstein constraints 



with transformation 



If a multi-skyrmion describes atomic nuclei upon quan- 
tization, 

it has to be quantized as a boson or as a fermion 
whether B is even or odd. This requirement is imple- 
mented in the form of Finkelstein-Rubinstcin (FR) con- 
straints The FR constraints for the rational map 
ansatz was constructed in 13] and [s^ and applied to 
predict the ground states of skyrmions up to B = 22. In 
this section, we shall apply the FR constraints for the 
rational map ansatz directly to our axially symmetric 
multi-skyrmions and obtain their ground states. 

Following the notation in ^33] , let 17 be a rotation by 
a around n followed by an isorotation by f3 around N. 
Then the FR constraints can be defined as 



exp(— zan • J) cxp(— i/3N • T)^/j = XFR{g)^ (49) 



where 



XFR{g) 



1 if contractible 
— 1 otherwise. 



and, J and I are space-fixed spin and isospin operators 
respectively. "0 is the wave function which transforms 
under a tensor product of rotations and isorotations. In 
particular, a closed loop is noncontractible for odd B and 
contractible for even B, which is consistent with spin 
statistics. Consequently, quantum numbers / and J are 
half-integers for odd B and integers for even B. 

In order to construct the ground states for a given 
baryon number B, let us define N{L(a, (3)) as a homo- 
topy invariant for a loop L generated by rotations by a 
and isorotations by f3. Then, for the axially symmetric 
rational map of degree i?, it is given by js^ 



N{Lia,P)) 



B 
2^ 



{Ba - P) . 



(50) 



It can be shown that N (mod 2) determines if the loop 
is contractible or not in the same sense as B (mod 2). 
Therefore, N (mod 2) gives the FR constraints for each 
generator of the symmetry group of the rational map. 

The axially symmetric rational map with degree B is 
given by 



R{z) = 



1 



(51) 



There are two symmetric generators for this rational 
map. One is a rotation by a followed by an isorota- 
tion hy P — Ba. Substituting it into (|50|) . one obtains 
N{L{a, Ba)) — 0. The FR constraints for this loop is 
thus given by 



iniL3-BK3)j, _ 



1p — 1p . 



(52) 



where we introduced the body-fixed spin (L) and isospin 
(K) operators related to the space-fixed operators by or- 
thogonal transformations. The other symmetry is C2 



, Riz) 



R{z) 



(53) 



This corresponds to a = j3 ~ tt and hence N{L{tt,7t)) 
B{B — l)/2. The FR constraints for this loop is 



(54) 



In the following we construct the ground states con- 
sistent with the derived FR constraints H52|l and (|54|l for 
_B = 2 — 5 with axial symmetry. 



• B = 2 

We find the FR constraints 

-i7r(L3-2K3),/, _ 



i/j = i/; 



,-i7r{Li+Ki) I _ 



1p = —Ip . 



(55) 
(56) 



This gives the ground state as |J, La)]/, X3) = 
|1,0)|0,0). 



• B = 3 

We find the FR constraints 



(57) 
(58) 



This gives the ground state as ||,|) Ij)^)- 

• B = 4 

We find the FR constraints 

This gives the ground state as |0,0) |0, 0). 



• B = 5 

We find the FR constraints 



(59) 
(60) 



-0 = ■)/) 



-i-iT(Li+Ki) I _ 



(61) 
(62) 



Tp = tp . 

This gives the ground state as | J) |/) 

|Z 5\ 11 1\ 
I 2 ' 2 / I 2 ' 2 /• 



Thus, for even _B, the axially symmetric solitons are pos- 
sible candidates of the ground states of B atomic nuclei 
as is the case of the deuteron and ^He while for odd B 
they emerge only as excited states. 
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IV. NUMERICAL ANALYSIS 
A. Eigen Equations 

In this subsection, we show the numerical analysis of 
the eigen equations in detail. To solve the eigenequation 
of the form, 

[—icy. • V + /3Af (cos F{p, z) + i"f^T ■ n sin F[p, z)](/)^(a;) 

^E^(j,^{x), (63) 

we introduce the deformed harmonic oscillator spinor ba- 
sis which was originally constructed by Gambhir et al. in 
the relativistic mean field theory for deformed nuclei ■ 
The upper and lower components of the Dirac spinors are 
expanded separately by the basis as 



tions, and ff„^ , are the associated Laguerre polyno- 
mials and the Hermite polynomials with the normaliza- 
tion constants 



inr + \uj\y. 



(68) 



These polynomials can be calculated by following recur- 
sion relations 



(64) 

where #a(35, s, r) arc the eigefunctions of a deformed har- 
monic oscillator potential 



Voscip, z) = ^Mujy + ^Mcjy , 



and defined by 



^aix,m,) = ^^H(p)0„^_(z)e-^x^, 



(65) 



(66) 



with 



^^:}ip) = A^tl(V^)l"le-^"''P^LH(app2) 

Tlr — 0, 1, 2, • • •, Nrmax 
max 

+ 1 or 0, 2, • • • , 2Af, 



and 



X4 



X- 



(67) 



depending on if the eigenvalues of the third components 
of the spin to^ (isospin to,-) takes +1 or —1. The func- 



x^Kix) = nL^{x) -{n + m)LZ^,{x) 



and 



Hn+i{x) - 2xHn{x) + 2nzHn-i{x) = 
—Hn{x) = 2nHn-i{x) 



(69) 
(70) 

(71) 
(72) 



where constants ap and can be expressed by the os- 
cillator frequencies as 



Mujo 



MoJz 



(73) 



h ' ^ h 

which are free parameters chosen optimally. The Nrmax 
and Nzmax are increased until convergence is attained. 
The parity transformation rule of <I>q is given by 

$„(p,(^ + ^,-z;s,<) = (-l)"+"-$,(p,<^,z;s,i) (74) 

where 

i?n.(-V^2) = (-l)"^i/„^(V^z), (75) 

has been used. The parity is + for uj + = odd, and — 
for (jj + Hz = even. 

There are two sets of the complete basis for each parity. 
One is the natural basis with = 0+, 1^, 2+, • • • , for 
odd B and = , | , • ■ ■ for even B. Another is 
the unnatural basis with = 0^, 1+, 2^, • • • , for odd B 

and = ^ ' f"*^' f ' ' ' ' ^'^^ even B The natural basis 
is given by 



q(0) a(l) 

«5If/3(o),M*/9(o)(a;, Ts) + ij^am-.i^^mi^' Is) 



/3(0) 



J 



"(2) 
\ /3(2) 



ix,U) + ^fai3).M^a{3)ix, is) \ 

J 



(2),M*a(2) 
a(2) a(3) 

«E5'/3(2),M*/3(2)(a;, Ts) + i^5'/3(3),M*/3(3)(a;, is) 



m 







where 



and 



a(0) 


= {rir, 


Uz 


odd. 


UJo 


= K3 


-1/2- 


TOw/2} 


m 


= {nr,nz 


even 


Wo 


= i^3 


-1/2- 


TOw/2} 


a(l) 


= {"-r, 


Hz 


even 


Ul 


= K3 


+ 1/2- 


TOw/2} 




= {nr,nz 


odd. 


UJi 


= K3 


+ 1/2- 


TOw/2} 


a{2) 


= {n-r, 


Uz 


even 


UJ2 


= K3 


- 1/2 + 


TOw/2} 


/3(2) 


= {nr,nz 


odd. 


UJ2 


= K3 


- 1/2 + 


TOw/2} 


«(3) 




nz 


odd. 


U3 


= Ki 


+ 1/2 + 


"T-w/2} 


/3(3) 


= {nr,nz 


even 


W3 


= K3 


+ 1/2 + 


TOw/2} 
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The unnatural basis 0^"'' is given by replacing, a <-!■ /3 in 



B. Matrix elements of the eigenequation 

By using the natural and unnatural basis, the eigen- 
value problem in Eq. can be reduced to a symmetric 
matrix diagonalization problem. 

Let us calculate the matrix elements of the Hamilto- 
nian below. For the kinetic term 



we have 



Ck • p 



cr • p 
cr • p 



{^c.(o)xi\M cos F{p,z)\^ 
= J pdpdzM COS F{p, z) 

= — j pdpdzM cos 0(/O, z) sa\F{p, z) 

x4r;'(p)0n.(^)</>L":'(p)0n'.(^) 



(79) 



(80) 



and 



we have 



(*a(0)|o- -PN^/S'CO)) 



1 



n'-l 



and 



(*a(0)|o- •PN0/3'(1)) 



xe 



V 9p p dtp/ 



(77) 



(wo > : = cjQ + 1 > 0) 

Tii-n' — 1 y 

(tJo < : = Wo + 1 < 0) 



(78) 



In the natural basis, quantum numbers {nz,n'^) 
takes values (1, 2), (3, 4), • • • for the upper part and 
(1, 0), (3, 2), • • • for the lower part. In the unnatural ba- 
sis, {nz,n'z) = (0, 1), (2, 3), • • • for the upper part and 
{riz, n'z) = (2, 1), (3, 2), • • • for the lower part. 
For the potential term 



f3M{cosF{p, z) + ij5T • iisin F{p, z)) 

cos F{p,z) iT ■ n sin F{p, z) 



= M 



-IT • nsh\F{p, z) — cos F{p,z) 



{^aiQ)xi\MiT • nsinF(p, z)\i^ ^2)Xd) 
= —J pdpdzM sinO{p, z) sin F(p,z) 

x0ir;i(p)0».(^)0y(p)</'n'.(^). (81) 



Other elements can be calculated in the same manner. 
In Appendix^ we shall present these calculations of the 
matrix elements in more detail. 



C. Matrix Elements for the Moments of Inertia 

To compute the matrix elements of the moments of in- 
ertia, we shall evaluate (n|Ti|m), {n\T3\m) and {n\Ji\m). 
For (n|Ti|m), only following elements survive 



{^i%,xi\r,\^il 



= hj c''2;0L?'(p)0n.(^)e-^"°^0Lf' (p)</'ni(^)e^"^'^ 
l(p)0„,(z)e— ^^0l:;^l(p)^„,(z)e-> 



2) 



in) 



l(p)0„,(z)e— ^^0l:;«l(p)^„,(z)e'-ov' 



{<%)Xi\ri\^': 



{u) 



(82) 
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For {n\T-i\m), 

(*a(0)Xik3|*a'(0)Xf.) = 1^3 | (1) X^) 

= (*/3(0)X^I'r3|*/3'(0)X^) = (*/3(l)XfjT3|*^'(i)xf,) 

(*a(2)Xdk3|*a'(2)Xd) = (*c(3) Xd ^3 I *a' (3) Xd) 
= (*/3(2)Xd|T3|*/3'(2)Xd> = (*/3(3)Xdk3|*/3'(3)Xd> 

= — f^n^n' 5„^„' Jjj-gjj-' . (83) 



For (n|Ji|TO), we shall evaluate (n|[_H'(L/'Q'^), Ji]|m) in- 
stead with the replacement of 144|l . 



(n|[i/({/o^^),Ji]|™) 

= (n| [/3A/(cos z) + i75T • n sin F{p, z)), Zi] |to) 
= — (n|Zi[/3Af (cosi^(p, z) + i75T • nsinF(p, z))]|to) 



where 



h = 
h = 



-—I 
2 



V 9z So odtjD/ \ 



d_ 

dp 



d 



+ 1- 



d 



dz dp pd^p' 



dz dp pdip) V dz dp p dip 



dip 



(84) 



are the components of the angular momentum operator I in cylindrical coordinates. Using eas. H84|l . one obtains 



( a a .z a \ ^, . 

x<|pi""i+i-''i+^4?'Kp')4f'(«pp') 

(nzi7„^_i(Va7z) - -iJ„^+i(^/arz)^i?„;(ya7z) 
+-H"«,(\/aIz)(^n^iJ„^_i(Varz) - iif„^+i(Varz)^ 

[((2n,. + l^ol - ay )Ll7l(ay ) _ 2(n. + |c.o|)4?Ii(«pP'))4f' («pP') 
+Llr;l(app^)((2< + |..^| - app2)4f I(app2) - 2(n; + |c.^|)l|:^^ ,(«pp2)) ± Lj^;! (apP^)4?' Kp') 

X Hn, {^/a^z)zHn'^ (v^z)| COS i^(p, 2)'5;f3zp 1 1 



(85) 
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3 z d \ 

z— =Fi-— sinF(p, z)cos8(p, z) 

oz op pOipJ 



K^/3'(0)/ 



wo| + |(^o|-l 



X (ya7z)ziJ„/^ (\/arz)| saiF{p, z) cos e(p, z)(5;^^^ i ± i 



(86) 



- 2^ ^*^^) sinF(p,z) sine(p,z)£ 



-trriip 



l*^/3'(2) 



(n^Hn^^l{y/Q^.z) - ^Hn^+l{y/c^z)jHn'^{y/a2z) 



-Hn^{^/a^z)[ji'^Hn^-i{^/a^z) - -Hn^+ii^/c^z) 



((2n. + Ic^ol - a,p^)Lt:Ko^pP') " 2K + ko|)4rii(«pP'))4f' («pP') 



x^^„,(V^^)-z-f^"L(\/^-2) r sinF(p, z) sine(p, 



(87) 

Other elements can be obtained in the same manner. In Appendix ^ we shall present these calculations in more 
detail. 



r 



D. Numerical Convergence 



In this subsection, we show the convergence of the soli- 
ton energy Estatic with respect to T) K3 with the basis 
number fixed, and U) the number of the basis (discretized 
momentum number) with K3 fixed. 

FigllUI shows the case of T) with B = 3. As can be 
seen, the energy is almost convergent at the {K3)max — 
10. In the case of U), the energy is not perfectly conver- 
gent up to Nrmax = A^zmaa; = 22 (scc FigHIl)- Therefore, 
all our results have 1 — 2 % uncertainty. 



CONCLUSIONS 



In this paper, we studied axially symmetric soliton so- 
lutions with S = 2 — 5 in the chiral quark soliton model. 
The one-quark spectral flow analysis indicates that the 
number of diving states from positive continuum to neg- 
ative coincide with the baryon number. As is shown in 
Table m the valence quark spectra contain double degen- 
eracy, realisng lower energy than non-degenerate states. 
Therefore, our solitons are stable although they are not 
necessarily minimal energy one, which confirms that they 
are good saddle point solutions. 

Upon quantization, we computed zero-mode rotational 
corrections to the classical energy. The study of the 
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3570 
3560 
3550 
3540 
> 3530 
3520 
3510 
3500 
3490 
3480 
3470 



TABLE V: -B = 4, mass spectrum up to i < 3, j < 5, fcs < 1. 



FIG. 10: Convergence of (Tfa)™ 

tion with Nrmax — Nzrnax — 10. 



for B = 3 calcula- 



Classification 


ihj, kz) 


Parity 


Mass[A/eF] 




(0,0,0) 


+ 


4753 


(0,2,0) 


+ 


4807 




(0,4, 1) 


— 


4808 




(1, 1,0) 


+ 


4813 


(2,0,0) 


+ 


4879 


AC S3) 


(1,3,0) 


+ 


4904 




(2,2,0) 


+ 


4934.2 


(0,4,0) 


+ 


4934.3 




(2,4, 1) 




4935 


N ^A^Ps) 


(1,5,1) 




4941 


A^Si) 


(3,1,0) 


+ 


5025 


N 'ACSs) 


(1,5,0) 


+ 


5046 


^A^Si) 


(2,4,0) 


+ 


5061 


''N ACSs) 


(3,3,0) 


+ 


5115 


N 3A(9P5) 


(3,5,1) 




5152 


N ^A^Ss) 


(3,5,0) 


+ 


5278 



3750 - ■ 
3700 - 
3650 - 



3500 - 




3450 I I I I 1 I I I 1 I I I 1 I I I I I I L. 

4 6 8 10 12 14 16 18 20 22 

Discrete momentum number 



FIG. 11: Convergence for the basis number - EataUc- 



{Kz)max is fixed by 6. 



TABLE IV: B — 2, mass spectrum up to i,j < 3, fcs < 1. 



Classification {i,j,kz) Parity Mass[MeV] 

NNCSi) (0,1,0) + 2264 

NNCSo) (1,0,0) + 2290 

7VA(=^P2) (1,2,1) - 2399 

NA(:'S2) (1,2,0) + 2477 

NACSi) (2,1,0) + 2528 

AA(^5'3) (0,3,0) + 2576 

AA(i5'o) (3,0,0) + 2730 

AA(^P3) (2,3,1) - 2762 



Finkelstein- Rubinstein constraints indicates that the ax- 
ially symmetric solution with even B has the same quan- 
tum number as the physically observed nuclei. These 
results are shown in Table HVI and IVI Some of the states 
may be observed in experiments. For odd B, the con- 
straint of 6*2 in Ea. H54l) seems to assure the validity 
of the ansatz. Indeed, it provides the ground state as 
i^j = 1/2 for S = 3 and as i = 1/2, j = 3/2 for 5 = 5, 
which exactly coincide with physical observations. This 
seems to make sense since in the minimal energy con- 
figurations with discrete symmetries, the solutions tend 
to have i = j = 1/2 due to their shell- like structure. 
However, unfortunately the constraint in Ea. H52l) forbit 
such states. Consequently, the axially symmetric solitons 
with odd B emerge only as excited states. The resultant 
lowest state is E = 3657 MeV with i = l/2,j = 5/2 for 
5 = 3, and is E = 6591 MeV with i = 1/2, j = 7/2 for 
5 = 5. 

Recently, we also studied classical multi-baryonic sys- 
tems with discrete symmetries in the CQSM and found 
larger degeneracy of the quark orbits than of the axi- 
ally symmetric j^, H^l- For example, triply degenerate 
bound spectrum is obtained in the 5 = 3 tetrahedral 
soliton background. The interesting point is that the cor- 
responding energy of the soliton is ~ 210 MeV which 
is higher than the axially symmetric, = 173 MeV (see 
Table QJ. Likewise, for the 5 = 4 minimal energy soliton 
with cubic symmetry, the valence quark spectrum shows 
four-fold degeneracy with E ~ 170 MeV while for the 
axially symmetric, E = 106,232 MeV. Thus, although 
the degeneracy of the spectum indicates the stability of 
the solutions, other factors should be also taken into ac- 
count in regard to minimization of their classical ener- 
gies. More detailed discussions on this subject will be 
made elsewhere. 
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APPENDIX A: EVALUATION OF THE MATRIX 
ELEMENTS 



In this appendix, we shall present detailed calculations 
of the matrix elements, / (Px(f>jj^H(j)^. 

For the kinetic term of the elements, one gets 
1 



2tt 



ar(^n'^i/„'^_i(Varz) - ii?„^+i(^arz) 



xTV, 



^nj.n'^^ujQUjQi ^ riz ^ n'^ 



2 "^"^'+1 



(Al) 



For (n,,<) = (1,2),(3,4),--- or (n,,0 = 
(0, 1), (2, 3), • • • , the righthand side of Ea. ljAlll becomes 



Nn' 



(A2) 



For (n„<) = (1,0),(3,2),--- or (n„ <) = 
(2,1),(4,3),---, 



2 -"rij.n'^'JuJoujQ jy- 



(A3) 
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Also, 



(*„(0)|cr.p|z*^,(i)) 



xL^:^\^a,p') - 2« + \u;[\)L^:^Uap')) . (A4) 



For ojQ > {oj'i = 070 + 1 > 0), the righthand side of 
Ea. (|A4|l becomes 



= Sn^n'^ J pdp{^f'^° + 'p'"^°e-''-P A^;r:Cr' 



1 



^0.0 + 12 



-Vrir + Ci'O + l^n,.n;, + \/n^(5„^ - 1„;, ) 



where 



+1 



iV"f V«r + C^O + 1 ' 



(A6) 



wo + 12 



= y^Tlr +UJO + l<5n^n' , (AT) 



(A8) 



have been used. 



Likewise, for uq < = ujq + 1 < 0), 

x(L,--Kp2)-L;,;_!,Kp2)) 



-1^0-1 



1 , 



2 '^nj-n'^ — 1 



(A9) 



where 



A^r°^< = V^^<^n.n-l, (AlO) 

A^r-A^^r"'!^ -v^^^^'5„,<' (All) 

have been used. In the same manner, elements, 

(A13) 

can be analytically calculated. 

Since cr-p does not contain isospin operators, following 
elements vanish 

(*a(0)k •PN*/3'(2)>(*a(0)k -PN^/S'IS)) 



(A14) 



For the potential term in the Hamiltonian given by 

^^1^ cosF(p,z) ir • hsmF{p,z)\ 
-ir • h sin F{p, z) — cos F{p,z) J ' 
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only following matrix elements survive 

(*a(0) 1^^ COS z) I *„,(o)) 

= j pdpdzM cos z) 

(*a(o) N^'^T • " sin z) I *^/(o) ) 
= — j pdpdzM cos Q{p, z) sin F{p,z) 

(*c<(o)NA/t •nsinF(p,z)|*^/(2)) 
= — J pdpdzM sin Q{p, z) sin F{p, z) 



Other elements of the potential term can be calculated 
in the same manner. 



(A15) 



(A16) 



(A17) 



APPENDIX B: EVALUATION OF THE 
MOMENTS OF INERTIA 



In this appendix, we shall give the detailed derivation to obtain matrix elements in Ea. H85|l . 



For the first term in Ea. (|85(l . one obtains 



^ J d^xcl,t°\p)cbnA^^)e-'^°''e^^^(^p-^^ cos F{p, 



= / pdpdz4r;l(p)^„Jz)(p^cosF(p,z))0l:;°l(p)0„,(z)<5;^3^i;^,±i 



(Bl) 



Performing partial integration with respect to z in Eq. (|B1 



pdpdz^l,^! (0„^_ (z)0„, (z)) cos F{p, z)5, 

- / pdpdzN\;:'j'\N^;;'°\^)\'^°\+\'^'°\p\'^"\+\^«\+^e-^ 



n^Hn^^iiy/a^z) - -Hn^+iiy/a^z))Hn>Jy/a^z) 



cosF{p,z)5k -riK'± 



(B2) 
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Similarly, performing partial integration with respect to p, one obtains 



-±-^z|-cosnp,z))|*a„p 



(o)i^ '"'["dp 



^ I d3a;0KI(p)</,„^(^)e--o0e±iv'(^^eosF(p,^))0i^°l(p)</,„,^ 

I pdpdz4>t'\p)'i>n. iz) (^^ COS Fip, Z)) </.Lf (P)-/-™; (^)5^3^ 1 ^± 1 
- J dpdz^ (pcAk?' ' (P)) <^„. (-2)^</'ni (^) COS F{p, Z)5j,^^ij,,±i 

X ((2n. + |a;o| - a,p')L\l^\{a,p'') - 2K + ko|)4?IiKp'))4f' Kp') 
+47l(app2)((2< + - a,p2)4f IKp^) - 2« + |a;^|)L|:;^,(a,p2)) 

XNn^Nn'^e'""''^ Hn,{^z)zHn'^{y/^z) cos F{p, z)6K^^iji'±i . 



Finally, we reach the final answer 



/ d d .z d \ s 
\.Pd-z-'d-p^%d^)''''^^P^'\ 



l^a'(O)/ 



(B3) 



.l-ol+l-ol-i [((2n, + |a;o| - a,p'')L\l"\{a,p') - 2K + ko|)4?IiKp'))4f' Kp') 



X Hn, {^/a;z)zH„,^ {^/a^z) | cos F{p, 



z)5_ 



(B4) 



